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Abstract 

The order parameters of dynamical chiral symmetry breaking in QCD, the dy- 
namical mass of quarks and the chiral condensates, are evaluated by numeri- 
cally solving the Non-Perturbative Renormalization Group (NPRG) equations. 
We employ an approximation scheme beyond "the ladder", that is, beyond the 
(improved) ladder Schwinger-Dyson equations. The chiral condensates are en- 
hanced compared with the ladder ones, which is phenomenologically favorable. 
The gauge dependence of the order parameters is fairly reduced in this scheme. 

1 Introduction 



The dynamical chiral symmetry breaking plays an important role in the field of particle 
physics. In particular, the hadron dynamics gives the most important example of the 
dynamical chiral symmetry breaking phenomena. At present the low energy phenomenol- 
ogy of pions and kaons are well understood by regarding them as the Nambu-Goldstone 
bosons of the dynamical chiral symmetry breaking S77(3)l X S77(3)r — > SU(3)y caused 
by the non-vanishing quark pair condensate (qq) 7^ in QCD. 

There have been many studies about the dynamical chiral symmetry breaking. The 
ladder Schwinger-Dyson (SD) equations in the Landau gauge have been mostly and ex- 
tensively used.|l|-|l(J In the strongly coupled QED, the chiral critical behaviors were 
explored and the anomalous dimension of operator ipip as well as the critical exponents 
of the phase transition, were obtained. @-@ The SD approach has been applied also to 
various models beyond the standard model. @, §|, §J The chiral order parameters for 
QCD: the dynamical mass of quarks, the chiral condensate and the decay constant of 7r 
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mesons /„- were also calculated by using the improved ladder approximation in the SD 
equations. [0, || They incorporate the running effects of the asymptotically free gauge 
coupling constant in the ladder self consistent equations, thus called the improved lad- 
der approximation, and offered good results, even in quantities. However it should be 
noted that the improved ladder approximation has no theoretical justification, and it 
has been just an artificial "model". Therefore there has been no way to improve fur- 
ther the "improved ladder" until the Non-Perturbative Renormalization Group (NPRG) 



11, O, 0, 14) analyses shed new light on it. [15, 16] Furthermore it has been shown that 



the ladder SD equations are plagued with the strong gauge parameter dependence. || 
Also it is difficult to proceed beyond the ladder approximation so as to overcome this 
unpleasant problem. fllPl 

The NPRG is also a powerful analytical technique for the study of non-perturbative 
phenomena in various fields of physics. |l7j The NPRG has many favorable properties: The 
NPRG equations can be formulated exactly and there are systematic methods of approx- 
imation. The exact NPRGs are given in the form of the nonlinear functional differential 
equations for the Wilsonian effective action. Therefore, it is necessary to approximate 
them for the practical calculations. We expand the Wilsonian effective action in powers 
of derivatives and truncate the series at a certain order. We often use the lowest order 
of this approximation, the so-called Local Potential Approximation (LPA).[|18| Although 
the LPA looks like a very crude approximation, it gives very good results indeed, for 
example, for the second order phase transition in O(N) scalar field theories, compared 
with other non-perturbative methods, 1/N-expansion and e-expansion. [^0|1 

In previous papers, J15], [16|, ^TJ we studied the chiral critical behaviors mainly in the 
Abelian gauge theory with strongly coupled massless fermions. We proposed a set of 
the gauge independent NPRG equations describing the chiral phase transition. Gauge 
independence is achieved by taking in the non-ladder type diagrams to the f3 function 
of the four-fermi operators, as well as the ladder type ones. The gauge independent 
values of the critical exponent and of the anomalous dimension of the mass operator 
ipip were obtained. When we restrict the (3 functions to the "ladder parts", the NPRG 
equations exactly reproduce the ladder SD results. It was analytically shown in the 
ladder approximation that not only the critical behaviors but also the order parameters 
are identical to the results of the SD equations. It is remarkable that this exact equivalence 
between the ladder part NPRG and the ladder SD holds also for the improved ladder SD 
with the running gauge coupling constant. Thus our NPRG method has given for the 
first time a definite physical meaning to the improved ladder SD, and now we know how 
to improve the improved ladder, which we will challenge in this article. 

In Ref. [^TJ we proposed a new scheme of the NPRG equations incorporating the com- 
posite operators so as to definitely evaluate the order parameters of the dynamical chiral 
symmetry breaking. This formulation was also applied to the non-Abelian gauge theo- 
ries by taking account of the asymptotically free running of the gauge coupling constant. 
However the practical calculations of the order parameters were demonstrated only in the 
ladder approximation. Therefore the results obtained there should depend on the gauge 
parameter just as the ladder SD equations suffer. Note here that our calculational method 
of incorporating the composite operators should be distinguished from other schemes of 
constructing the effective meson theory at some scale in QCD, although they should be 
compared with each other. Also non-perturbative renormalization group analyses of QCD 
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with effective meson components are done in a scope of the hadronic matter. p2 



In this paper we evaluate the chiral order parameters of QCD in a new approximation 
scheme which is a minimal extension including the non-ladder type diagrams so as to 
compensate for the serious gauge dependence of the ladder parts. The results obtained in 
this scheme, however, are not completely free of the gauge parameter dependence. There 
are still a few sources in our approximation causing the gauge dependence. We examine 
the amount of the gauge dependence appeared in the order parameters. It is found that 
the gauge dependence in the observable quantity, the quark condensate (V>V)> i s fairly 
reduced compared with the results in the ladder approximation. 

The outline of this paper is as follows. Section [| is a brief review of the formalism 
upon which our work is based, the method of NPRG and its approximation. We present 
our model in section 0, where we show how to treat the infrared divergences occurring 
in the dynamical chiral symmetry breaking. In section |] we consider the origin of the 
gauge dependence in the ladder approximation, and construct a new beyond the ladder 
approximation which should reduce the gauge dependence. In section || we describe 
the practical calculation procedures and we present numerical results. We discuss the 
remaining gauge dependence of the order parameters in our approximation in section [f| 
Summary will be given in section [7]. 



2 NPRG Equation and its Approximation 



There have been known several formulations of the NPRG. ITT], O, O, [TJJ In this paper we 



take the Wegner- Houghton (W-H) equation. [12] First let us briefly review its formulation. 



The starting point is the Euclidean path integral with the controlled momentum cutoff 
A(t) = e-*A : 

f A(t) 

Z = J P0exp[-5 cff (0;t)], (1) 

where S'eff is called the Wilsonian effective action. The NPRG equation describes how the 
Wilsonian effective action S e s should change as the higher momentum degrees of freedom 
are integrated out. It is obtained by reducing A(t) infinitesimally with fixing the partition 
function Z. Simultaneously we rescale the momentum variables and the fields by cutoff 
A(t), since the change of the dimensionless quantities are of our physical interest. We 
obtain the following differential equation, 

dS eS [(j>;t] r d D p ( & \ 5 

- DS cS - / — -— - <pi{p) Dj, - p^— — — 



dt Ln J\p\<i (2ir) D V dp»J WiiP) 

1 f d D p f 6S C{{ ( 6 2 S cS V 1 aSeff 



2dtJ {2-K)D\8(j )i {jp)\8 ( j> i {jp)84> j {-p)J 6(f>j(-p) 



5 2 S, 



- strl HwJf (2) 

where D is the space time dimension, is the dimension of (f) including its anomalous 
dimension, the second primed integral denotes integration over the infinitesimal shell 
modes of momenta e~ dt < p < 1, and the prime in the derivative indicates that it does 
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not act on the 5 function in -^-S c s. The subscript i represents every Lorentz and internal 
symmetry indices. This equation is known as a sharp cutoff version of the NPRG, and is 
called the Wegner-Houghton (W-H) equation. |TJj It is inevitable to approximate them 



for the practical calculations. We expand the Wilsonian effective action in powers of 
derivatives. We employ the Local Potential Approximation which is regarded as the 
lowest order of this derivative expansion. Any derivative couplings are dropped except 
for the fixed kinetic terms, 

S eS = /^{^,(p)0, + K ff (0)}, (3) 

where Kij(p) is a matrix of the canonical kinetic terms, and V e s is called the Wilsonian 
effective potential. As a simple example, we consider a theory of one scalar (p and one 
Dirac fermion if) and its conjugate if). The matrix K(p) in (0, if), if>)-spsice is written as 

/ p 2 \ 
K(p) = -if . (4) 

V o -%i> o J 

In this approximation the W-H equation is reduced to a nonlinear partial differential 
equation for the Wilsonian effective potential V e ^{4>, *)>P 



= DV„ - + i f'^str In U s + , (5) 

dt * dfa 2 J (2tt) d V d<\>id<\>i) 

where D'^ denotes the canonical dimension of field 0. The anomalous dimension of field 
vanishes because the kinetic term is not renormalized in the LPA. We should notice that 
any newly generated operators with derivatives are ignored in this approximation. We 
take account of the generated operators which do not depend on the external momenta. 
Definitely speaking, we evaluate amplitudes of such local non-derivative operators by 
setting all external momenta to vanish. 

This partial differential equation may be solved numerically. However, actually it is not 
easy to get its solution with enough precision. Besides it would not be practical for more 
complicated models. Therefore here we expand the effective potential into polynomials in 
field 0. By this approximation we solve a system of coupled ordinary differential equations 
for the various coupling constants. 

Let us characterize these approximations from the viewpoint of the NPRG formulation. 
The basic logic of approximation in the NPRG formalism is to restrict the theory space 
to a subspace of the original full theory space. Namely the approximation in the NPRG 
formalism is to analyze the NPRG equation projected onto a subspace, which is actually 
finite dimensional so as to get results numerically within finite computation time. In 
order to improve the approximation, we enlarge the subspace, step by step, expecting the 
results will converge to certain values. The above two approximations, the local potential 
approximation and further the polynomial expansion, are just two consecutive steps of 
the subspace projection. It should be noted that the ladder approximation itself can not 
be regarded as projection to any subspace, and therefore it has some pathological features 
indeed. 

Finally we should mention an intrinsic problem of the NPRG formulation. The mo- 
mentum space cutoff is indispensable for almost any formulation of the NPRG. It comes 
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out that the NPRG does not manifestly respect the gauge invariance. There have been 



several approaches for restoration of the gauge invariance. [23, 24, 25] Our purpose here 



is not to overcome this gauge invariance problem but to improve the improved ladder ap- 
proximation. Therefore we take a simple approximation scheme for the gauge interactions. 
We evaluate the j3 function of the gauge coupling constant by the one-loop perturbative (3 
function. Of course the running of the gauge coupling constant is automatically derived by 
the original NPRG equation. When we adopt a sub-theory space with lower dimenaional 
operators, then the NPRG equation effectively reproduces the perturbative renormaliza- 
tion group equation.^, |T5| Therefore as the fist stage, we suppose this level of the small 
sub-theory space, and we adopt the running gauge coupling constant controlled by the 
perturbative j3 function. We should notice that any newly generated operators including 
the gauge fields are irrelevant in this scheme. Therefore the Faddeev-Popov ghosts are 
also irrelevant. 



3 NPRG equations for the dynamical chiral symme- 
try breaking 

Now we apply Eq. (0) to QCD with three massless quarks. We take the local potential 
effective action, 



s cS [^, v, A; <] = J d x [ v es(^, ^ t) + v# - g4)i> 

(6) 



where a is the gauge parameter, and ip denotes massless triplet quarks. Here, as mentioned 
at the end of the previous section, the gauge coupling constant g is supposed to follow the 
one-loop RG equation. We start with the general form of the effective potential consistent 
with the chiral symmetry S[/(3)l x S77(3)r and the parity. Let us first consider four- 
fermi operators. We regard the operators corresponding to the Fierz transformation as the 
identical operators. Furthermore we shall not consider the flavor and/or color changing 
multi-fermi operators. Then there are two independent four-fermi operators in our theory 
space. 



= (^) 2 + (^ 75 ^) 2 = -i{(^) 2 -(^757^) 2 }, 
2 = $7mV0 2 + $757^) 2 - (7) 



To approximate the NPRG equation, we specify a subspace in the full theory space. 
Here we take a subspace spanned by polynomials in scalar operator 0\ only up to some 
maximum power nmax, 

nmax /~i f±\ _ 
n n 

The NPRG equation for the scalar four-fermi operator is obtained from the diagrams 
in Fig. g, 

^_ 2G2 + _L (G2)2 + g G2 + i£. (9) 
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Generally the NPRG equation is a set of coupled equations of various operator. In this 
case, however, the NPRG f3 function of four-fermi operator 0\ consists of 0\ itself and 
the gauge coupling constant g. Therefore the NPRG equation for 0\ operators can be 
solved without recourse to other higher multi-fermi operators. When ignoring the running 
of the gauge coupling constant, there is a fixed point for the flow of (D\ operator given by 
the zero of the right-handed side of Eq. which is nothing but the critical point of the 
dynamical chiral symmetry breaking, and we have two-phase structure of the standard 
ferromagnet type phase transition. [|15|, [16| The strong coupling phase is the symmetry 
breaking phase, where the four-fermi coupling constant diverges at a finite scale. 

We make the gauge coupling constant run according to the asymptotically free f3 
function. The flow diagram in the G2 — g plane in QCD is depicted in Fig. |[ There 




>0< + T + 









Figure 1: The NPRG function of the 
four-fermi operators. The diagrams in the 
solid box correspond to the ladder part, 
and the diagrams in the dashed box are 
crucial for the gauge independence. Possi- 
ble other diagrams are ignored in our ap- 
proximation. 




Gauge Coupling Constant 



Figure 2: The flow diagram in the G 2 -g 
plane in QCD. The phase structure does 
not appear, and all flows diverge at a cer- 
tain finite scale. 



appears no phase boundary. All the flows diverge at a certain finite scale. Also we see 
the renormalized trajectory in the flow diagram which assures that the bare four-fermi 
interactions are irrelevant to the infrared physics. This behavior suggests that the entire 
region is supposed to be in the broken phase of the chiral symmetry. This is due to the 
infrared slavery behavior of the QCD gauge coupling constant and it is believed to hold 
naturally. 

Renormalization group flows have been also analyzed by the SD equation method pfifl. 
where fixing the quark mass obtained, the relation between the four-fermi coupling con- 
stant and the gauge coupling constant is calculated assuming the cutoff dependence of 
the gauge coupling constant. This procedure is justified within the SD formalism and 
it actually gave something resembling to the results in Fig. |2|. However there are criti- 
cal differences between these two calculations. The main difference comes from the fact 
that in the NPRG formalism the bare four-fermi interactions are turned out to be irrele- 
vant, while in the SD formalism there is no mechanism to automatically generate effective 
four-fermi interactions by the gluon exchanges. 

Correspondence between the divergence of the four-fermi operator and the dynamical 
chiral symmetry breaking is not trivial. All the coupling constants in the polynomial 
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expansion keep growing in the infrared and diverge at some finite scale. This corresponds 
to the fact that at this scale the Wilsonian effective potential exhibits a non-analytical 
behavior at its origin, which is actually observed as a jump of the first derivative by 
direct analysis of the Wilsonian effective potential using the partial differential equation. 
This singularity at the origin clearly shows the spontaneous symmetry breakdown since 
it assures the non- vanishing magnetization at zero external field limit. 

This singular behavior invalidates the renormalization group calculation of the evo- 
lution of the effective potential with polynomial expansion at the origin. Introducing a 
composite operator corresponding to the order parameter (V>V) enables us to carry out 



the calculation of the Wilsonian effective potential even in the infrared region. Our 
theory space is extended to include the composite operator 0. First we introduce the 
composite operator as an auxiliary field in the original path integral without changing 
the dynamics. The partition function Z in this extended theory space is written as 



VipVipVAexp 
VijjVijjVAVcpexp 



d x C OTg [ip, ip, A 



d A x\C mg [4>^,A} + ^{(i ) -y^) 2 



(10) 



Here we abbreviated the counter part ('ipi^ip) 2 for simplicity and it should be considered 
that the original chiral symmetry is kept actually. Then the bare lagrangian is modified 
as follows: 

£ = C OTg [ij,ij,A; t] + I(0-^) 2 



\f; v f; v + ^(<V^) 2 + </# - 9 4- + \<P 2 + ji 



Even in case of omitting the chiral symmetric representation, this lagrangian has the 
discrete chiral symmetry: 

$ -> 7 5 ^, V -$7s, -» -0- (12) 
Then we may analyze the NPRG equations for the following Wilsonian effective potential: 

As = -F^ + ^-g^ + V^a; t), (13) 

Kff(0,<x; t) = G o (0; t)+ £ -G„(0; t) a n 

= G o (0; t) + G 1 (4 >] t)a + ^G 2 (0; t)a 2 + ■ • • , (14) 

where the notation a = x/np is introduced. Note again that we are actually working with 
the chiral symmetric effective potential and we take a particular direction of the scalar 
operator a condensation to get the effective potential represented by Eq. flT4]). Since we 
do not consider the propagation of composite operator modes, Eq. ( |I4] ) is enough to define 
the NPRG evolution of the chiral symmetric system. 

In this formalism, it was shown that the chiral condensate (ipip) is proportional to 
the minimum position of the scalar potential Gq(4>), denoted by (0), which is the order 
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parameter of the dynamical chiral symmetry breaking. 

<^> = -(0), (15) 

y 

Then the dynamical mass of quarks S(0) is given by 

=G 1 (0=(0)) (16) 

(T = 

In the usual argument of introducing the auxiliary field, the four-fermi interaction is re- 
moved away from the action by tuning the Yukawa coupling constant y. Then this action 
may be regarded as the gauged Yukawa system with the compositeness condition. || There 
are several studies of the dynamical chiral symmetry breaking by using this realization. 
However, we should notice that our purpose of introducing the auxiliary field is abso- 
lutely not to eliminate the four-fermi interaction. Also our results obtained by Eqs.(|l5D 
and fljjD are turned out to be independent of the Yukawa coupling constant jy,pl] since 
it should not change the dynamics at all. 



E(0) 



dV f 



eflf 



da 



4 Compensation of the gauge dependence 



As noted before, the ladder part NPRG exactly reproduces the results obtained by the 
ladder SD equations. Namely the results by the ladder part NPRG depend on the gauge 
parameter a strongly, like the ladder SD approaches. In order to improve the gauge de- 
pendence, we must proceed beyond the ladder approximation. So it is required to develop 
a non-ladder extended approximation in the course of the systematic approximation of 
NPRG. 

To begin with, we discuss the origin of the gauge dependence in the ladder approxi- 
mation. Let us consider a set of diagrams summed up by the ladder approximation. For 
that purpose we define the "massive" quark propagator, 



1 



m 



(17) 



where 



7710, a 



QT/ C) nmax /~t / a nmax 

da da ^ n ^ 

Gi(<j>) + G 2 {<P)a + G 3 (0)a 2 + G 4 (0)cr 3 + • • • 



Using the Feynman diagrams the "massive" quark propagator is represented as shown in 
Fig. [3[ The ladder part f3 functions are defined by summing up a set of diagrams which 
do not contain any crossed ladder type diagrams (Fig. f|). 

We now consider the Abelian Ward-Takahashi (WT) identities assuring the gauge 
independence. To satisfy WT identities we must sum over the diagrams for the S-matrix 
at any given order. When the gauge boson is inserted at a certain point along the fermion 
line, we must sum over all possible insertion points. Let us consider the simple case of WT 
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(—)'=(—)'- (_ + JL + i^ + ...) 

1 I 



Figure 3: The "massive" quark propagator. The deep full line in the left-handed side is 
a "massive" fermion propagator. The pale full lines are massless fermion operators and the 
dashed lines are the auxiliary fields 4>. 




Figure 4: The ladder part (3 function. These diagrams do not contain any crossed ladder 
type diagrams. The wavy lines are gluons, and the deep full lines are the "massive" quarks, 
and the pale full lines are external quark operators. 



identity involving two fermions and two gauge bosons in g 2 order. The sum of diagrams 
in Fig. H is gauge independent at on-shell. The right diagram of Fig. j|, called the crossed 
diagram, is not involved in the ladder approximation. This is one of the reason for the 
strong gauge dependence in the ladder approximation. Actually as seen in Fig. [I], adding 
the crossed box diagram to the (3 function has wiped out the gauge dependence of the 
critical behaviors. [[15] 

Now we consider to generalize the crossed box diagrams in Fig. [I] for the f3 functions 
of the higher multi-fermi operators. First we define a corrected vertex, 

" W^ 7 57/3 + m<T , (19) 



p2 _|_ m 2 



which is composed of "two" diagrams, ladder and crossed, using the "massive" quark 
propagator (Fig. |6[). Therefore the corrected vertex itself comprises an infinite number of 
diagrams. Then we replace double vertices in Fig. § with the corrected vertex and sum 
up the diagrams just as the ladder part (Fig. [7]). 




lip' |4p vcy 

Figure 5: The set of diagrams satisfying the Ward-Takahashi identity which involves two 
fermions and two gauge bosons. The wavy lines are gauge bosons and the full lines arc fermions. 
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~& - W + "X" 

Figure 6: The corrected vertex. The wavy lines are gluons, and the deep full lines are 
"massive" quarks, and the pale full lines are external quark operators. The curved arrows 
denote the direction of the shell-mode momentum p. 

= O + d + h 



+ M M + 

\ 

Figure 7: The j3 function in our new approximation. The wavy lines are gluons, and the deep 
full lines are "massive" quarks, and the pale full lines are external quark operators. 

Just as in Fig. f|, every diagram in Fig. |7] contributes to an infinite number of coupling 
constants in the polynomial expansion. Accordingly, a certain diagram with n external 
quarks contributes to the (3 functions of coupling constants with more than n external 
quarks. For instance, the (3 function of mass G\ (the generalized Yukawa coupling) is 
given by the first and the second diagrams in Fig. |7]. Also the f3 function of the four-fermi 
operator comes from the first, the second and the third diagrams, etc. Note that the third 
diagram has a symmetry factor of 1/2. 

Now we discuss the way of "projection" of the newly generated operators onto our 
target subspace. Here we consider the effective potential V e g composed of the polynomials 
only in the scalar operator 0\. Therefore we pick up only the parts proportional to (l) n 
in the spinor space from the generated total operators. For example, let us consider g 8 
terms in the eight-fermi beta function (Fig. |8|). The form of the generated operators are 
in general, 

£ c a u^r>i) x $ 2 rfy 2 ) x (^ 3 r c ^) x (^r^), (20) 

a,b,c,d 

where T's are the 16 independent matrices for spinor indices and C a bcd are amplitudes. 
We should consider all the parts proportional to (ipli/j) 4: from the generated operators. 
Instead of picking up all of them, however, we take only one simple combination of T's, 

(Vil^i) x faith) x OMV-a) x (toO ■ (21) 
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Figure 8: The newly generated g 8 operator. The numbers denote the suffixes of the quark 
helds (see Eq. ©). 

This part is represented in the left figure in Fig. [| Of course there are another combi- 
nations of r's contributing to (tplip) 4 term, for example, contraction shown in the right 
figure in Fig. || which is omitted. We adopt this approximation since the left figure in 
Fig. |9| exactly coincides with the ladder part approximation, when the corrected vertices 
are replaced with the ladder type. We now challenge the minimal extension of the ladder 
approximation. We should mention that this restriction of the way of picking up a fields 
does not correspond to a general procedure of projection in the systematic approximation 
method of NPRG, and it might cause a problem. 




Figure 9: The way of picking up a. The left figure corresponds to the ladder part approxi- 
mation, when the corrected vertices are replaced with the ladder type. We ignore the way of 
picking up a like the right figure. 

Now let us consider the non-Abelian effects. The quarks belong to the 3 dimensional 
representation. For example, g A term in the four-fermi (5 function from the ladder type 
diagram (Fig. 10) is identical to the Abelian estimate, except that we must add the color 



Casimir eigenvalue, 



b 




Tr (T a T h T h T a ) = 3C 2 (S) : 



Fig. [T(J. The gauge group factor of the box diagram in the four-fermi /3 function. 

As for the crossed ladder type diagrams (Fig. [11]), there appears the commutator term 
if abc in addition to the Casimir term. Here we simply ignore the commutator term if abc 
and take account of the C2(3) 2 term only. 
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T a c f 

/ \y \ ^ ^Tr(T a T b T a T b ) 

\ A / 3C 2 (3) 2 + if abc Ti(T a T b T c ). 

jib j,a 

Fig. 0. The gauge group factor of the crossed box diagram in the four-fermi /3 function. 

Of course this is a violent truncation which might break the gauge independence and it will 
be discussed later. Due to this additional approximation the non-Abelian nature is ab- 
sorbed into the Casimir factor which defines the effective gauge coupling constant C 2 (3)g. 
With this effective gauge coupling constant, the f3 function of the effective potential in 
QCD is evaluated just as in the Abelian case. 

Now we are able to write down general formulae of g 2n terms in the NPRG equations 
as follows: 




,2n 



8tt 2 2"" 2 n! (1 + m 2 ) n 1 + 5 n>1 + 8 n , 2 



x 



hi 



"•(3 + ar )+ g(-)> (a) , (B _ a) , 



(2 + A k )m 



n-2k 



a 



For example, the j3 function for the four-fermi operator reads 



dG 2 



2G 2 4 
9 



G\ 



2vr 2 l + Gj 
2 3 + a 



2G\ 
1 + Gl 



4tt 2 l + G\ 



Go 



2G\G 2 
l + Gj 



+ 



16tt 2 



G\ 



6 - (3 + a 2 )G 



Let us compare above equation with the ladder approximated one, 



dG 2 



-2G 2 + 
9 



G\ 



2vr 2 l + Gj 
2 3 + a 

Go — 



1 - 



2G\ 



4tt 2 1 + Gl 



2G\G 



l + Gi 

2 



G\ 



9 4 (3 + a) 



32tt 2 



1 + G 2 



1-Gl 



(22) 



(23) 



(24) 



The gauge dependent terms in Eq. (^3|) should be almost compensated by the gauge 
dependent anomalous dimension of the quark fields when we go beyond the LPA, which 
will be discussed later. 
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5 Numerical Calculation and Results 



Now we describe how to get the chiral order parameters in QCD with our approximation 
scheme. We work with the Wilsonian effective potential defined in Eq. ( |i4|) with some 



finite highest powers nmax, and we numerically integrate its NPRG equation. The NPRG 
equation is defined by the j3 function given in Eq. fl2"2]), that is, we take only the quantum 
loops of the quarks and gluons and not of the scalar composites. The initial effective 
potential is taken from Eq. (0). During evolution the scalar field is fixed to be a 
certain value just as an external source field. The gauge coupling constant is set to follow 
the one-loop perturbative (3 function with three flavor quarks. We take the QCD scale 
parameter Aqcd to be 490 MeV and we also adopt the same infrared cutoff scheme of 
the gauge coupling constant divergence as in Ref. ||, since our results should be first 
compared with the previous ladder SD results. 

Integrating the NPRG equation, the effective potential finally stops to move except 
for the canonical scaling behaviors, where the cutoff scale has been lowered well below the 
quark mass scale. Then we get the scalar potential Gq(4>) at the fixed value. To solve 
the NPRG equation scanning the fixed <ft value, we obtain the scalar potential function 
Go(<f>) and its minimum point {(f)). Then we estimate the chiral condensates and the 
quark mass using Eqs. (|T5|) and flUf). The chiral condensates obtained above should be 
regarded as the bare operator condensation at the initial highest cutoff scale. It should 
be renormalized through the standard procedure to show the renormalized condensates 
at 1 GeV scale. 

We show the results (a = case) in Fig. 12. First of all we check the nmax dependence 
of the results. Though there still remain some small fluctuations, we may claim that we 
have obtained the results of our total subspace of nmax = oo. We have checked also that 
the dependence of the initial cutoff A should be negligible, that is, our results are assured 
to be on the renormalized trajectory and we can regard them as those of the infinite initial 
cutoff limit. 
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Figure 12: The chiral condensates and the dynamical mass of quarks with their truncation 
dependence in the Landau gauge. The dashed line is Aqcd- Non-ladder results are enhanced 
compared to the ladder results. 



Now we compare our non-ladder extended results with the ladder ones. The ladder 
results exactly coincide with those of the ladder SD equation, which assures the total 
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consistency of our calculational machinery. The chiral condensates and the quark mass are 
both enhanced by including the non-ladder contributions. Though we do not argue here 
in detail about the phenomenological implications of this enhancement, this enhancement 
is actually favorable for phenomenology since our setting of the QCD scale parameter has 
been recognized to be much higher than the current estimate even considering that it is 
the one-loop (3 estimates. 

The gauge parameter dependence of the results are depicted in Fig. |T^. Compared 
with the ladder results, the improvement of gauge dependence of (ipip) is clearly seen 
in the left figure of Fig. |T3|. On the other hand the right figure of Fig. [T3| shows that 
the gauge dependence of S(0) still remains a lot even in the non-ladder approximation. 
We understand the different situations between these quantities as follows. The chiral 
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Figure 13: The gauge parameter dependence of the chiral condensates and the dynamical 
mass of quarks are plotted for cases a — (Landau gauge), a = 1 (Feynman gauge) and a = 2. 



condensate (V>V) 1S a measurable physical quantity. And it should not depend on the 
gauge. On the other hand the dynamical "mass" of quark X(0) is an off-shell quantity 
and is not directly related to a measurable quantity; therefore, it may depend on the 
gauge. It should be noted also that the quark mass E(0) strongly depends on the infrared 
cutoff scheme of the gauge coupling constant divergence while the chiral condensates do 
not. 



6 Issues of the gauge dependence 

In this section we discuss the origin of the gauge dependence in the NPRG method. First 
of all, the Wilsonian effective action S e fi itself depends on the gauge, because it is not 
directly related to any measurable quantities. Therefore the NPRG equations (or the (3 
function) describing the evolution of the Wilsonian effective action also depend on the 
gauge. Furthermore even at the infrared limit, the effective action is not totally gauge 
independent except for the on-shell quantities. For example the effective potential is not 
gauge independent except for the position of the minimum. Thus in general the (3 function 
which depends on the gauge parameter finally gives the gauge independent results only 
for physical quantities at the infrared limit. Actually in our approximation scheme of 
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evaluating the effective potential, there is no way of erasing all the gauge parameter 
dependence in the (3 function. 

We discuss here the gauge dependence due to the approximation we adopted, that is, 
the Local Potential Approximation. It ignores any corrections to the derivative couplings 
including the kinetic terms, and therefore no anomalous dimension is taken into account. 
This seems to be the largest source of the gauge dependence. We will report elsewhere 
the results taking account of the quark anomalous dimension, where we will see the 
reduced gauge dependence. Before getting these new results, we may evaluate the physical 
quantities as follows. In the one-loop approximation the quark anomalous dimension is 
proportional to a, and therefore we have vanishing anomalous dimension in the Landau 
gauge a = 0. Therefore we may claim that the Landau gauge results in our scheme are 
most significant and they would be very near to the coming results with quark anomalous 
dimension. Then our main results should read, 



1/3 



Aqcd 



= 0.512 ±0.014, (25) 

non— ladder 



which is compared with the previous ladder results 

(^) 1/3 



Aqcd 



= 0.439. (26) 

ladder 



There are other subtleties of the gauge dependence due to the LPA. In the above 
calculation we have done further approximation to ignore some parts contributing to the 



(3 functions. In Eq. (22) we omitted the non-Abelian commutator parts of the gauge 
effective vertex. This omission itself does not bring the gauge dependence. Rather it 
"hides" the gauge dependence of the LPA. Consider the four-fermi amplitude, for example. 
Such commutator parts should be summed up with diagrams in Fig. |14] (second derivative 
part of the gluon field, d^F^ip^ip as an operator form) to generate the gauge independent 
(on-shell) four-fermi amplitudes. This situation is quite the same as the Penguin diagram 
to give the local four-fermi effective operators. Of course we also have to add all related 




Figure 14: The "Penguin" diagrams contributing to the four-fermi (3 function. 

diagrams in Fig. |l| and ghost diagrams to get totally gauge independent results with 
the properly renormalized gauge coupling constant. Therefore we have to take account 
of the derivative couplings d^F^ip^ip to compensate the gauge dependence appearing in 
the four-fermi box diagrams. 

All these extension requires the higher order derivative couplings in our sub-theory 
space. Then we have to proceed to use smooth cutoff scheme NPRG equations since the 
sharp cutoff NPRG equations suffer singularities when applied to the derivative couplings. 
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Figure 15: The diagrams containing the wave function renormalization of the gauge field in 
the four-fermi amplitudes. 

7 Summary and Discussion 

In this article we challenge a beyond the ladder calculation of the dynamical chiral symme- 
try breaking in QCD by using the non-ladder extension in the Non-Perturbative Renor- 
malization Group method. The ladder approximation of the NPRG Local Potential j3 
function has been integrated to give exactly the same results of the (improved) ladder 
Schwinger-Dyson equation for the chiral condensates and the dynamical mass of quark 
S(0). Extension beyond the ladder has strong motivation of reducing the inevitable gauge 
dependence of the ladder approximation. 

We add non-ladder diagrams to the NPRG ladder (3 function, trying to reduce the 
gauge dependence of the physical results. We develop a set of (3 functions using the 
effective gluon vertex defined by the sum of the ladder and the crossed couplings. We nu- 
merically solve this new f3 function to get the chiral condensates and quark mass function 
at zero momentum. They are enhanced compared with the previous ladder results, which 
are favorable phenomenologically. Also we evaluate the gauge parameter dependence of 
our results and find it is fairly reduced compared to the ladder case. 

We stress here again that our results are the first results in the long history of ana- 
lyzing the dynamical chiral symmetry breaking in gauge theories, which goes beyond the 
(improved) ladder equipping with a systematic approximation method. This is realized by 
quite a new viewpoint of the NPRG method for the dynamical chiral symmetry breaking. 
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